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Abstract 

In relativistic heavy ion collisions, the cross section of heavy hadron 
production via the combination of a heavy quark with a light one from 
the quark-gluon matter can be factorized. It is the convolution of twist-4 
combination matrix elements, the parameters corresponding to the parton 
distributions of the quark-gluon matter, as well as the hard partonic cross 
section of heavy quark production calculable in PQCD. These parton dis- 
tributions and combination matrix elements are functions of a scaling vari- 
able which is the momentum fraction of the heavy quark w.r.t. the heavy 
hadron. In the same factorization framework, the combination matrix el- 
ements appear in other 'simpler' processes and can be extracted. Taking 
them as inputs, comparing with data from RHIC and (future) LHC, we can 
get the parton distributions of the quark-gluon matter just as the similar 
way we get those of nucleon, pion, or photon, etc.. 



1 Introduction 

One of the main purposes of (ultra-)relativistic heavy ion collision experiments, 
e.g., those in RHIC at BNL and LHC at CERN, is to produce and study a de- 
confined phase of quarks and gluons, usually referred to as Quark-Gluon Plasma 
(QGP), under the extreme conditions at high temperature and/or high density. 
QGP is an important prediction of Quantum Chromodynamics (QCD) and is sup- 
posed been existing in the early universe shortly after the Big Bang [1]. In the 
past few years, the hot dense quark-gluon matter has been produced at RHIC, 
while experimentalists are still working on more signals [2]. At such a stage, 
aware of the particular complexity of the experiments, one can be more conscious 
than ever that to get more knowledge of QCD at extreme conditions and/or to 
'simulate' the state of the universe at early time on collider, one should go further 
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than just to 'discover' QGP. That is, to carefully and quantitatively measure 
it. In fact, it is widely agreed among experimentalists and theorists that more 
detailed measurements are needed before drawing a definite conclusion on what 
is produced at RHIC [2, 3]. There have been many marvelous probes suggested, 
which have played the main role in the study of the hot dense quark-gluon matter 
at RHIC [2, 4]. On the other hand, new methods with unambiguous quantitative 
relation between the final state particles and the partons of QGP (or other states 
of quark-gluon matter produced in the collision) are still worthy to be looked 
for. In this note, we study the possibility of measuring the Quark-Gluon Matter 
(QGM) via the heavy hadron produced by a heavy quark (charm or bottom) 
created in the hard collision and combined with the partons from the QGM. The 
data from RHIC themselves are the suggestion for this idea. 

Recently, some 'unexpected' phenomena on hadron production in Au-Au col- 
lision at RHIC ([2] and refs. therein), e.g., different p± suppression between 
mesons and baryons and the L v 2 scaling', challenge the 'energy loss + fragmenta- 
tion' picture of hadron production and are well explained by various combination 
models [5, 6, 7]. Besides, these models provide a clue to measure the partons 
in the QGP/QGM because the production rates of the hadrons are proportional 
(via convolution) to the distribution function of the quarks (gluons x ) in the 
QGP/QGM. Hence one can derive, or at least infer the 'parton distributions' 
from final hadrons. The work of [6] is an example. The authors extract the 
quark distribution from the pion data within the framework of their model, and 
predict or explain the production of the others. 

The quark (re) combination models are of a long history. They can date back 
to three decades ago. From then on, different kinds of quark (re) combination 
models are presented in application to hadron production in various high energy 
processes[8]. Therefore it is not surprising that several different combination 
models are suggested to explain the RHIC data as cited above and new works 
are emerging [9]. Unfortunately, such a condition puts forward the problems and 
uncertainties: Both the distributions and combination rules of these models are 
different from each other, so one is difficult to gain the unambiguous knowledge 
on the content of the QGM. In other words, the 'parton distribution' inferred in 
this way is model-dependent. If one wants to 'measure' the parton distribution of 
QGP/QGM in the combination production processes, one has to fix the unique 
'combination rules'. This seems beyond approach for light hadrons now. 

However, the combination of a heavy quark (charm or bottom) with light ones 
from the QGM can provide more clear information, and is free of model depen- 
dence. The main points are: 1) The combination production can be factorized. 
The cross section formula (e.g., that of the inclusive process A + B — > Mq + X, 
here A, B, Mq denote the colliding nuclei and the produced heavy meson, re- 

1 Thc phenomenological models cited above only consider the valence quarks of a meson or 
baryon now. 
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spectively) is the convolution of the hard sub-cross section of the heavy quark 
production (A + B — > QQ), the combination matrix elements and the parameters 
corresponding to the parton distributions of the QGM. The combination matrix 
element and the parton distribution of the QGM are expectation values of field 
operators on certain particle states, so they are model-independent and process- 
independent. 2) The spectra of heavy quarks can be calculated by PQCD (now 
full NLO corrections are available and tested in Tevatron [10].). 3) Within the 
same factorization scheme, the combination matrix elements which describe the 
probability of a heavy quark and a light one of specified momenta to form a heavy 
hadron, appear in other more 'simple and clean' processes, such as e + e~ anni- 
hilation, DIS and hadronic collisions, so can be extracted for the experiments. 
From the above points, it is clear that in the cross section formula, taking Point 
3) as inputs, comparing with heavy ion collision data, we can get the parton 
distribution of QGM. 

A simple comparison of the 'combination probe' with the hard probe of the 
energy loss/jet quenching can help to understand the idea in this note. We 
do not expect to find advantages or disadvantages, but we find that they are 
complimentary for each other. Let's first see the electroweak probe (photon, 
W + ~ , Z°) on nucleons, in which case, the single photon (or other gauge bosons) 
approximation is enough. Imagining that if on the contrary, we have to sum 
the electroweak coupling to all orders, i.e., multi- interaction between the electron 
(neutrino) and the proton by exchanging multiple gauge bosons is important, we 
lose the unambiguous definition of W^ u . This is the case for a hard jet in QGM, 
where multiple soft strong interactions between the hard parton and the QGM 
have to be taken into account. So, it is hardly possible for the hard jet to play the 
same role of the electron (neutrino). However, for the combination process of a 
heavy quark with the light one from the QGM, the cross section is proportional to 
the parton distribution function of the QGM. It is more attractive in the results 
of this note that, for certain Q, which is the CMS energy of the partonic process 
of heavy quark production, the heavy meson of momentum K just comes from the 
heavy quark with momentum fraction z = 2 , hence probes the light quark in 
QGM with momentum k — (1 — z)K. This is very similar to the Bjorken scaling 
in DIS. 

Heavy quarks are mostly produced from the initial parton scattering (in this 
note we do not discuss the possibility of creation of heavy quark pairs from the 
QCD vacuum fluctuation at very high temperature), they evolve in space-time 
and interact with the QGM. The measurement of jet quenching can give impor- 
tant information of the space-time structure of the QGM, e.g., the correlation 

2 This is exact for the 2^2 partonic process. For final states with more partons, the relation 
is kept as long as we take Q as the invariant mass of the heavy quark pair. As we shall see 
in Section 2, this relation comes from the on shell condition of the 'freely-fragmenting' quark 
and the 4-momentum conservation for the quark pair system, which are also the origin of the 
Bjorken scaling. 
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length or the space scale of the medium, L [11]. On the other hand, the formulae 
in this note are the result of the integration on the whole space-time, so they are 
only sensitive to the local phase space (momentum) distributions of the partons. 
In fact, from the discussions in the following sections, we can see that the in- 
duced radiation can be included into the evolution of the partons or the matrix 
elements ( l z scaling violation'). So the energy loss and combination are both 
sub-processes in the heavy hadron production, and both are necessary in predic- 
tion of the experimental data. This note only concentrate on the combination 
mechanism. 

The outline of this note is: Section 2 introduces the derivation/factorization 
of the cross section, taking the case of heavy meson as an example . In section 3 
we show that the combination matrix elements appear in and can be extracted 
from other simpler processes. Section 4 gives the numerical estimation of the 
combination process cross section with the input of assumed values of the quark 
distribution of 'QGP' and the recombination matrix elements (since they are still 
beyond available from experiments now), only to demonstrate the feasibility of 
the formulae. Section 5 is for conclusions and discussions. The combination of 
heavy quark with gluons and heavy baryon production are left for later works. 

2 Combination of heavy quark with light one 
from QGM 

In the following of this note, we take A+B — > D+X as an example. Here D refers 
to any anti-charm meson. We only consider the contribution of c + q to D, other 
possible combination processes such as c + g to D are assumed negligible. The 
X includes the associated produced c quark and all the other particles from the 
nucleus-nucleus A, B interaction. In this section, we derive the inclusive invariant 

dN AB ~ 

differential production number 2E—^- of D. The subscribe C of N denotes the 
combination process. (E, K) is the 4-momentum of D. The production of charm 
mesons can be treated in the same way as the anti-charm mesons. 

To describe the light quark from QGM, we should find ways to represent 
the 'external particle source'. Similar as the works on energy loss [11, 12], we 
employ an external field to describe the interaction with QGM. We can see in the 
following, this external field, together with quark field operators, appears in the 
matrix element corresponding to the quark distribution in the QGM. Same as [13], 
we choose the external field a vector V^ xt proportional to n M = (0, l,0j_). This 
can be understood as gauge fixing and is easy to factorize the Dirac indices. We 
would emphasize here that the external field can present hot dense quark gluon 
matter or QGP, as well as cold quark matter (nuclei), as in [13]. If one calculates 
the matrix elements by, e.g., lattice QCD, one should identify the concrete form 
of the external field and discuss the details of the differences between hot and 
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cold QGM. However, here we only use it in the formal expression of the matrix 
elements, which is to be probed in experiments. So V ext is just a 'symbol' here. 

From the above discussions, the interaction Hamiltonian for quark and gluon 
fields is extended as: 

i?/ = #(4+y)# (i) 

Here A is the normal gluon field and V is the external field. The strong 
coupling constant g s is absorbed into the gauge fields. In this note, wherever we 
write the gauge field obviously, we always adopt this convention. One can easily 
see that this is the same interaction Hamiltonian in describing the energy loss in 
QGM by induced radiation, which has been discussed comprehensively in previous 
works, e.g., [11, 12]. In this note, we only concentrate on the combination. The 
jet quenching processes in our frame work can be taken as radiation corrections 
(i.e., to consider higher orders of the perturbative expansion of the S-matrix) in 
the QGM environment and will be treated as the evolution of the matrix elements 
(z scaling violation) in later works. 

If we assume that the distribution functions of the intrinsic heavy flavours in 
the initial nuclei are vanishing, the lowest order contribution for the recombina- 
tion process comes from O(gf) in the perturbative expansion of S-matrix: 

S (3) = ^pCf J ' d i x l d A x 2 d 4 x 3 T^(x l )A(x l )m(x l ) 

x y(x 2 )4(x 2 )y(x 2 )*(x 3 )Y(x 3 )y(x 3 ), (2) 

where summation on colour and flavour indices as well as the indices in spinor 
space are indicated. 

Now we take the annihilation partonic process qq — > cc as an example to 
illustrate the derivation, while the total result can be obtained by summing all 
kinds of partonic processes. Let the corresponding terms of the Wick expansion 
act on the initial nuclear state \AB > (We will discuss the Glauber geometrical 
formulae [14] in the end of this section.) and final state \DX >, employing the 
space-time translation invariance, we can isolate the 5 function corresponding 
to the total energy- momentum conservation and get the T- matrix element. The 
cross section is, then 



x < AB\i,{x A )T c ^i,{x^{x 3 ) y(x 3 )i/j(x 3 )V(x 5 )T c ^ 4 V(x 5 )\DX > 

x < DX\^(x 2 )T c ^^(x 2 )i}(0) y^^ix^T^^x^AB > . (3) 

In the above equation, T c is one half of the Gell-Mann Matrix, 4F represents 
the incident flux factor and discrete quantum number average. Summation on 
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repeated indices is indicated. The capital \l/ is for the heavy quark field and the 
lower case for light quark fields. 

The following thing is to factorize the cross section in the framework of 
collinear factorization. The cross section can be written as 

(4na s )*C f dS 1 
a = —[F—IW)i W 7 ^ () 

with 

= J d 4 x 4 e~ iqX4 < A^ixi)^ 1 (0)|4 >< B\4>^(x 4 )4> ai (0)\B > 

x T^^ + CA^S). (5) 

This is the same W^ v for Drell-Yan process, which gives the distribution of initial 
partons. We have employed the translation invariance and integrated on xi, which 
gives q — q' . 

On the other hand, 

D = f d \ K f d 4 x 2 d 4 x 3 d 4 x 5 e- lk ^e lk ^Y 

x (T^' + mh.W < 0|(^(x 3 ) y{x,))%\X h >< X h \(M(0)m)Z\° > 
x < ()r-(.r :i )vl^ (,,) D >< 5|*j,(x 2 )^(0)|0 > . (6) 

To get the above expression, we have written the total final state produced 
in the A B collision as \DX >= \DXcXh >, where the X represents all the 
particles except the c c produced by the hard interaction and Xh for the assemble 
of particles produced by all the interactions except the above hard one. We have 
the corresponding field acting on the c quark final state with momentum k', 
k 5 = q — k'. The summation on X has been eliminated by the completeness 
condition. 

In (4), C is the colour factor of the partonic diagram except the external leg 
combined into D. The colour part will be clarified following. 
The W^ u can be conventionally written as 

x < B\tr(^i,(y)m)\B > tr(^^-Y), (7) 

which is the intended factorized form. Some of the above variables are: Ai = 
P 1 f y_,A 2 = P 2 ~y+, V = (y + ,U~,0±), Q 2 = q 2 = x x x 2 s, s is the CMS energy for 
the nucleon-nucleon system whose partons collide and produce the heavy quark 
pair. 

The factorization for D^ u is more complicated. It can be written as: 
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/ 



cF>K 

: " ; dzK + dz l K + {27i)5{k' 2 -m 2 ){2^f5\h + k- c -K) 



(2ir) 3 2E 



x trdnifi + m)-f u —)\ k , =g - zK +Y^ 
A x h 



X 



x < 0|ir(^ r V jo (^)^(0))|5 >< 5|tr(^^(x 2 )^°(0))|0 > (8) 

In the above equation, the colour indices in the partonic final states and the 
distribution functions are summed. So the colour indices in the combination ma- 
trix elements should be averaged (|), which is similar as the case in fragmentation 
function. We do not separate the colour-singlet or the colour-octet contribution 
in the matrix elements, but sum them together. The reason is that the partonic 
cross section is the same for the colour indices j belonging to 1 or 8 states (since 
the other parton comes from an un-correlated source), and that for the parton 
distribution in QGM, it should be the same whether j belong to 1 or 8 states. 
We have taken the external field proportional to n M , which select only the "+" 
component since fi = fii£± jh — 0. 

In the part corresponding to the partonic sub-process, for the sake of factor- 
ization, we have done the collinear expansion for the momentum of the c along 
the "+" component of the momentum of D. Here one notices that the coordinate 
system is different from that for initial states W^. The z direction is along the 
the momentum of the anti-charm meson. 

In Equation (8), we write the combination matrix element (Row 4, 5) formally 
to be analogous to that in other processes (see Section 3). z and zi seem not 
restricted to be z + z\ = 1. However, the 5 functions in the first row sets the 
restriction. At the same time, the integral in the combination matrix elements 
/ d A ki acts on the matrix element corresponding to the quark distribution in the 
QGM (Row 3) as well as the 5 function in the first row. These show that we have 
not finished the factorization. To get the factorized form, we notice 



J d ^5\k l + k- c -K) 
I 



dK ~K K± 5 + (kf + kt-K+)5-5l 

— -2£5 3 (k c - + k ; -K) 

2E V c ' J 2K- 



-5+{kt + kt-K+). 



(9) 
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Let the 3-dimension 5 function absorbed into the combination matrix element, 
we get the 'restricted' matrix element or the dimensionless combination function: 



x < QM^4> jo (x 5 )V(0))\D >< D\tr(^-V(x 2 )^(0))\0 > 
x 2E5 3 (k 5 + k; — K). (10) 

The integral of z\ in the first row of Equation (8) have given z + z\ = 1. 
/ dz5(k' 2 — m 2 ) gives important result: z = 3 which is analogous to the 

Bjorken scaling variable. This means that for certain partonic CMS energy Q, 
the heavy meson with momentum K just comes from the heavy quark with 
momentum fraction z combined with light quark with momentum k\ = (1 — z)K, 
hence only probe this light quark in QGM. Such a conclusion does not depend on 
the special forms of the derivation in this note. In fact, the relation z = -^- q is 
set by the on shell condition of the heavy quark associatively produced with the 
one to be combined into the final state heavy hadron. Just like the DIS process, 
this is a physical condition which should be respected by any special forms of 
derivation. 

The cross section section now can be written as 



2E S§ = 4^E/^i^ 2 /a(^i)2/|(x 2 ) 

ab 



|2 1 (27T) 2 



X \ M a>r-^y^l)F(z, Zl )\ z+Zl=1 (11) 

Here f% and /g are parton distributions in Nuclei and should be treated 
by the Glauber geometrical formulae [14] . M is the mass of the heavy meson. 
\M\ 2 refers to the invariant amplitude square including all the coupling constant 
and colour factors for the partonic process ab —> c + x (where the momenta of 
external legs are modified and x to be considered as one particle). For example, 
for qq — > c + x, to the lowest order, \A4\ 2 is 

(^ as ) 2 cMY^Y Y) 7 trM ^' + m)lv ^ ] ' (12) 



3 This relation is exact for Q and K to be infinite while z fixed. It is a good approximation 
when the anti-quark to be combined into the heavy meson is on mass shell in the partonic 
processes. In this case, the four momentum of the anti-quark is (zK + , 2 ™ K+ , 0j_). In the heavy 

meson rest frame, it is easy to get z — ^ m +mV — jj- Go back to the initial parton 

CMS, use this approximation for the on-shcll condition (q— k) 2 = m 2 , we can get the relation. 
So, it is not just an approximation by taking M = m c = in K ■ q. 
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C is the colour factor. Though the quark mass term is vanishing, we keep it 
to show the origin of the formula. 

P(z{) can be understood as the distribution function of the parton (probed 
by the heavy quark) in the external source. P{zj) is also dimensionless: 



d A x 3 e-^\ k + =ZiK+ J2< X h \tr(mV(0)\0 >< 0\V(xMx 3 )^)\X h > 
x h 

(13) 

This can be understood as the expectation value on the state representing an 
assemble of particles produced in the A B collision denoted by \X h >. 

For numerical calculation, we give more discussions on the cross section for- 
mula. The heavy quark production sub-process in nucleus-nucleus collision is 
hard interaction and can be calculated in the binary approximation. The ini- 
tial state for this sub-process in Equation (3) then can be understood as nucleus 
modified nucleon. We use the Glauber geometrical formulae [14] to describe the 
distribution of nucleon in the nucleus, i.e., for a certain impact parameter b, the 
production (interaction) number Af is 

In the above equation, T AB {b) is the overlap function for nucleus A and B 

[14]. £ R is differential cross section for nucleon nucleon interaction calculated 
by Equation (11) with the initial state \AB >^ \NN >. dR is phase space 
element for the final state. The £ ~' on N indicates that the nucleon is modified 
in nucleus. Now In Equation (11), the parton distribution function is for partons 
in modified nucleon, The incident flux factor in Equation (3) is for nucleons and 
can be absorbed into the partonic cross section same as in the case of Drell- 
Yan process, with the above trace term (12) slightly modified to be the exactly 
partonic invariant amplitude square. 

Now the production number of D in the combination process is 



2E ^i = TAB{b) ^ l ^^A a ^0/2(^)^^^P(^)^(^^)|, +2i =l- 

(15) 

To get results for our relevance, we can integrate over b in central region. In the 
above equation, dX is the dimensionless invariant phase space for the '2-body' 
partonic final state c + x where x treated as one particle. This formula is also 
correct for higher order partonic cross sections. 
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3 The universal (process-independent) combi- 
nation matrix elements 



From the above section, The cross section of the D is dependent on both the 
combination matrix elements and the distribution of the light quark in the QGM. 
If we compare with data to extract the light quark distribution P(zi), one of the 
key inputs is the combination function F(z, zi), which is not calculable by PQCD 
and we should find ways to extract from more 'simple' experiments. This require 
that the cross section of a more simple process can be factorized and includes 
this parameter. The following is an example. Let's see the factorization and the 
complexity. 

It has been pointed out that, in hadronic interaction, the asymmetry of D 
meson in forward direction can be explained by the the combination of the initial 
parton with the charm quark produced in the hard interaction. Such a leading 
particle effect has been studied in [15], [16], both in the approximation m c — > 
oo. In such an approximation, the light quark has vanishing momentum, hence, 
qualitatively, the momentum of the D meson is approximately that of the charm 
quark, so that not possible to probe the momentum of the light quark. On the 
other hand, in [16], the authors also tried to give the combination matrix elements 
in the framework of collinear factorization, which is the same framework used in 
this note. The combination matrix elements there depend on 3 variables z±, 
z 2 , z 3 , seem not corresponding to the momentum fraction of the valence partons. 
However, starting from Equation (4) in [16], by taking into account the space-time 
transition invariance, we get the combination matrix elements with 2 variables 
corresponding to the momentum fraction of the charm and the light quarks, which 
is like those in the above section: 



r rl 4 h rl 4 h, k+ k + r 



(2vr) 4 (2tt) 4 



x <m k (0)^Qi(xi))\H Q ><H Q \Q i (0)^q j (x 2 ))\0>, (16) 



and 



J (2tt) 4 (2vr) 4 v K +! v K+ J J 
x <0\q k (0)^Q l (x 1 ))\H Q ><H Q \Q i (0)^ T q j (x 2 ))\0>. (17) 

The two parts of the combination matrix element, i.e., the double- vector part 
and the double-pseudo- vector part, should be separated here since the partonic 
cross sections corresponding to these two parts could be different. For the case 
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that the quark and the anti-quark from different sources respectively, as in the 
process in Section 2, these two parts can be put together and only the vector part 
needs consideration. 

The complexity lies in that, Equations (16, 17) are different from that in 
Section 2 by 2ES 3 (k s + k/ — K), without the restriction z + zi — 1. The reason 
is that in the process in [16], the light quark and the heavy quark can undergo 
hard interactions and in principle are not restricted on mass shell. We can also 
understand this from a different way. Equations (16) and (17) look like the 
combined distribution of two valence quarks in the heavy meson. In fact, rewriting 
the < 0|(---)i|#Q >< #q|(---)2|0 > to the form < H Q \(- ■ -) 2 \0 >< 0|(---)i|#q >, 
integrating the S functions and the exponential functions, taking |0 >< 0| = 1 
in the vacuum saturation approximation, we will get the form of the product of 
two parton distribution functions, each similar to that defined by Collins and 
Soper [17]. Then in a parton model at high energy, we can not require the sum 
of two parton momentum fractions equals one. Hence to get the inputs needed, 
we should find ways to relate the 'restricted' in Section 2 and the 'unrestricted' 
ones here. 

If the above matrix elements have been extracted from experiments, to get 
the restricted combination function, we start from, by denoting the Combination 
Matrix Element in Equations (16,17) as CME: 



CME = j^2E>5\K + K-K')y^2E> c 5\K-K) 

x ^E[5\^ - k^CME 
r d 3 K' ~ 

= J—F^z^z + zri (18) 

In principle, we should solve the integral equation and use the value of F on 
z + zi = 1 (K 7 = K) as our inputs. On the other hand, if the real world is more 
simple — as most models assume, F(z, z%\ z + z\) peaks around z + z% = 1, i.e., 
two valence quarks on mass shell with k c + ki ~ K — we can, approximately, fit 
CME as, 



CME - /w^n ^Tfa^ ^^) 

(«-0) = J w 2^ 3 (K-K')M 2 /(^,^) (19) 

So in the extreme/ideal condition, we can just have F(z, z t ) z+Zl= i oc ■ 
That is, because the distribution of F is a narrow peak, we use the average value 
of it in a reasonably small integral region of ^§r- 
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Equations (16, 17) will be applied in other processes and discussed in else- 
where, so that we can have more experiments to extract the combination matrix 
elements. In principle, when we accumulate enough number of data, especially 
from more than one process, the integral equation (18) can be solved. We just 
mention that we have discussed the combination process preliminarily in e + e~ 
annihilation, where a light quark 'fragments' into a heavy meson by combination 
[18]. 



4 Numerical estimation 

Since there are yet no data on the combination matrix elements and the heavy 
flavour particle production in heavy ion collisions, only to get a practical view 
of the discussions in this note, we give a numerical estimation of the transverse 
momentum distribution of the open charm mesons produced via combination by 
assuming the forms of the combination matrix elements and the parton distribu- 
tions of the 'QGP'. 

We make a simple assumption that the combination matrix elements are not 
sensitive to z and z\ and just take them as constant around z + z\ — 1. At 
the same time, because of no data, we can not specify any identified D meson, 
whatever D°, D ± or D* . 

There are many discussions about the possible distributions of the partons 
in the hot dense quark-gluon matter produced in RHIC. Especially, in all the 
combination models cited above, this is one of the key inputs. It is generally 
adopted that the distribution function of the thermal partons are exponential 
while those produced from the initial hard interactions are of minus powers. In 
the following we use some minus power, exponential and Gaussian distributions 
to see their differences. 

Because the absolute values of the parameters are unknown, as a demonstra- 
tion of the combination effects, we study the ratio of the D spectrum w.r.t. the 
charm quark spectrum, ^jp/^j and normalize the total ratio (integrated on pj_) 
to 1. 

Besides, there still some other things to be stated: 

1) The initial parton distribution function of nucleon in nucleus. Now there 
are several groups of nPDF's [19]. In this note, we just use the EKS one available 
in CERN PDF Library (version8.04) [20]. 

2) For simplicity, we use the LO partonic cross section. It is clear that to 
get the partial cross section of combination process from data, the contribution 
by heavy quark fragmentation should be calculated and subtracted from data. 
Such a work needs the NLO partonic cross section (see following discussion on 
the RHIC d Au data) and is in progress. 

3) After the heavy quark created, it inevitably go through the QGM and loses 
energy by radiation. From Equation (15), we can see that P{z{) ® F(z, z{) play 
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Figure 1: The ratio of the spectrum of D meson to that of charm quark, with total 
ratio normalized to 1. At the most left: the highest curve is for that the parton 
distribution of QGM, P(zi), is exponential; the lowest for P{zj) is Gaussian; the 
ones between them for P(z/)'s are of various powers. 

the role of Fragmentation Function (FF) of the heavy quark in this process and 
the radiation of gluons can be treated as the evolution of this FF (similar as the 
FF of a quark into the photon, see [21]). This leads to that, FF is originally 
function of z, now also of the scale fi 2 . In other words, the assumptions on P(zi) 
and F(z, z{) should be understood as the values on /z 2 = Q 2 . 
Figures (1) and (2) give the numerical results. 

From the most simple observation, because the heavy meson comes from the 
heavy quark combined with a light one, the momentum of the meson are always 
larger than that of the heavy quark. Such a fact is shown in the Figures that all 
the spectra are harder than that of the c quark. The curves are sensitive to the 
distributions, so that we can measure the parton distribution of the QGM from 
the heavy meson spectrum. 

Experiments are more interested on central rapidity region. We also give the 
ratio in the region \y\ < 0.5. The curves have the same property of those in the 
whole rapidity region. 

The RHIC data of D mesons in d Au collision show that their spectrum is as 
hard as that of the charm quark calculated at NLO. This may be an indication 
of combination (see [22]). Of course the external source for the light quark in 
this case is cold quark matter. Our calculation is consistent with these data 
qualitatively. The reason that we can not give any quantitative prediction or 
explanation is that the combination matrix elements are yet beyond approach. 



13 




-7 ' 
10 r 



2 4 6 8 10 12 14 16 18 20 



Figure 2: Same as Figure (1), except that the rapidity of the particles are re- 
stricted to central region (\y\ < 0.5). 

5 Conclusions and discussions 

The discovery of the extreme state of matter such as QGP is not the end of the 
story. We should continue to look for more ways to measure it in details. The 
combination process of heavy hadron production discussed in this note may shed 
light on such a purpose. 

The 'combination rules', suggested in various forms in different combination 
models, are defined by universal (process-independent) matrix elements in this 
note. They are model independent and can be extracted from other processes 
than the complex relativistic heavy ion collision. In the same framework, the 
'parton distribution function' of QGM (if allowed to be given such a 'standard' 
name), here also has an operator definition and is model independent. 

From the cross section formulae in this note, we can see that the matrix 
elements have the l z scaling behaviour'. We also have argued that the scaling 
violation can be explored by taking into account the radiation correction. 

The result of this note is just the beginning of a series of systematic works to be 
done, in theory as well as in experiment. This is just like the case in studying the 
parton distribution functions and fragmentation functions of hadrons. Among the 
works, the extraction of the combination matrix elements from various processes 
is basic and crucial for understanding the forthcoming experimental data. 

One thing many experimentalists and theorists are concerned is the collective 
behaviour of heavy quarks. The combination formulae above provide a baseline 
for observing this issue. If some collective behaviour is measured on final state 
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heavy hadrons, and if it only comes from the light quark which the heavy quark 
combined with, the value of the parameter, e.g., the v-i of the heavy meson should 
be \ or | times of light meson or light baryon, respectively, which is the value for 
a single quark. Only if the experimental value is larger than that, we can expect 
that the heavy quarks also have collective behaviour. This is a very interesting 
topic in study of the vacuum structure. In the T = vacuum, the heavy quarks 
can hardly be produced via soft interactions, e.g., 'tunneling effects'. However, if 
the temperature is extremely high and the QCD vacuum is changed, the heavy 
quark pair may be created and their collective effects may reflect the properties 
of the bulk. 

The author thanks Prof. X.-N. Wang for suggestion of studying this topic, and 
helpful discussions. He also thanks members of the Theoretical Particle Physics 
Group of Shandong University for encouraging and helpful discussions, especially 
Prof. Z.-G. Si for the discussion on combination matrix elements for the leading 
particle effects. This work is supported in part by the National Natural Science 
Foundation of China (NSFC) under Grant 10205009. 
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